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Long-range quasi-static gauge-boson interactions lead to anomalous (non-Fermi- 
liquid) behavior of the specific heat in the low-temperature limit of an electron or 
quark gas with a leading TlnT~^ term. We obtain perturbative results beyond 
CO [ the leading log approximation and find that dynamical screening gives rise to a 

I low-temperature series involving also anomalous fractional powers j'(3+2n)/3_ 

CSI ' determine their coefficients in perturbation theory up to and including order T^/^ 

Q^' and compare with exact numerical results obtained in the large-Nf limit of QED 

^ '. and QCD. 

^ '■ PACS numbers: ll.lO.Wx, 12.38.Mh, 71.45. Gm, 11.15.Pg 

^ ■ It has been established long ago fH] in the context of a nonrelativistic electron gas that the 

0> '. only weakly screened low-frequency transverse gauge-boson interactions lead to a qualitative 
^ I deviation from Fermi liquid behavior. A particular consequence of this is the appearance 
• of an anomalous contribution to the low-temperature limit of entropy and specific heat 
O . proportional to aTlnT"^ [US 01, but it was argued that the effect would be probably too 
2 ! small for experimental detection. 

I More recently, it has been realized that analogous non-Fermi-liquid behavior in ultra- 

p • degenerate QCD is of central importance to the magnitude of the gap in color supercon- 
p • ductivity 0, IE 0| ; cind it has been pointed out that the anomalous contributions to the 
D ! low-temperature specific heat may be of interest in astrophysical systems such as neutron or 
I protoneutron stars, if they involve a normal (non-superconducting) degenerate quark matter 

■ component. 

^ . So far only the coefficient of the aTlnT"^ term in the specific heat has been determined 

(with Ref. [3] correcting the result of Ref. by a factor of 4), but not the complete argument 
of the leading logarithm. While the existence of the TlnT~^ term implies that there is a 
temperature range where the entropy or the specific heat exceeds the ideal-gas value, without 
knowledge of the constants "under the log" it is impossible to give numerical values for the 
required temperatures. 

Furthermore, a quantitative understanding of these anomalous contributions is also of 
interest with regard to the recent progress made in high-order perturbative calculations of 
the pressure (free energy) of QCD at nonzero temperature and chemical potential 0, where 
it has been found that dimensional reduction techniques work remarkably well except for a 
narrow strip in the T-/i-plane around the T = line. 

In the present Letter we report the results of a calculation of the low-temperature entropy 
and specific heat for ultradegenerate QED and QCD which goes beyond the leading log 
approximation. Besides completing the leading logarithm, we find that ioi T / ^ <^ g <^ 1, 
where g is either the strong or the electromagnetic coupling constant, the higher terms of 
the low-temperature series involve also anomalous fractional powers T*^'^"'"^")/^, and we give 
their coefficients through order T'^/'^. 

Our starting point is an expression for the thermodynamic potential of QED and QCD 
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which becomes exact in the hmit of large flavor number Nf |9|, llO| , and which at finite Nf 
has an error of order g'^ in the regime T//i -C (?, 

P = NNf I + + > 
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where N = 3, Ng = 8 for QCD, and both equal to one for QED. Dt and Dl are the spatially 
transverse and longitudinal gauge boson propagators at finite temperature T and (electron 
or quark) chemical potential /x obtained by Dyson-resumming one-loop fermion loops, and 
Dvac is the corresponding quantity at zero temperature and chemical potential. 
Nonanalytic terms in the low-temperature expansion arise from the contribution 

f rdqo^^^^^^^j,^, (2) 
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The bosonic distribution function Ub = l/[exp(go/r) — 1] restricts the go integration to 
qo ^ T and T is assumed to be the smallest mass scale in the problem. Consistently 
dropping contributions proportional to in the pressure (T^ in the entropy), which for 
T / n <^ g are beyond our perturbative accuracy, it turns out that we only need the T — 
limit of the inverse propagator D^^, and only the lowest orders in qo/q and qo/ [i: 



ReD^^ = q^(l + 0{g 
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ImZ^^^ = {q' - ql) {l2^^' + + ^o^) ^^(2/. - q) (4) 

where = g'^Nj for QED and g'^Nf/2 for QCD. 

Keeping only the leading terms in the limit q^ ^ gives 
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Inserting this approximation into Eq. Q leads to the integral 
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Performing the go integration then gives the following contribution to the entropy 5* 
dP/dT (per unit volume): 



+0{T^/^). (7) 

While this reproduces the coefficient of the anomalous TlnT~^ term reported in j^, 
the coefficient under the logarithm as well as the suppressed 0(T^/^)-contribution are still 
incomplete. 

To complete the term linear in T, one has to perform an exactly analogous calculation of 
the longitudinal contribution, which involves 

and when inserted into ((H) contributes 

Finally, the remaining parts of (0), which do not involve the bosonic distribution function, 
yield 

Snon-ni, 9eS 2r 
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The latter contribution matches exactly the one from the standard perturbative result 
to order g"^ [11], while the contributions ((Tj) and ^ depend on having T/n <^ g. In this 
region, all of the contributions listed so far are negligible compared to the zero-temperature 
contribution ~ g'^^'^ in the pressure (which is only partially included in (^). However, by 
considering instead the entropy (and further below the specific heat), the above contributions 
become the dominant ones. 

Adding them all together, we obtain 

S-S, _ gl^^-T 4^ _ 2 + _ ^.{2) 
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+ c,/sT'/'' + cj/sT'/' + 0{T'), (11) 

where 5*0 is the ideal-gas value of the entropy per unit volume. 

To also obtain completely the coefficients of the terms in the low-temperature expansion 
which involve fractional powers of T we need to include more terms of and (@)) than those 
kept in (jni). A lengthy calculation, whose details will be discussed elsewhere, gives 

_ 8 22/3r(|)c(|) /3 

^5/3 9^3^11/3 ^9eSfi) , (12) 



_ 30 2V3r(f)C(f) 




FIG. 1: Transverse nfe-contribution to the interaction part of the low-temperature entropy density 
in the large- A'^y limit for the three values g^f^ = 1,4,9. The heavy dots give the exact numerical 
results; the full, dashed, and dash-dotted lines correspond to our perturbative result up to and 
including the TlnT"^, T^/^, and T^/^ contributions. 

Setting T//i ~ gl^^ with 6 > 0, one finds that the terms in the expansion (fTT|) correspond 
to the orders gl^^ \n{c/ ges), 9ea''^^^^^ ^ g^g^'^^^^^, respectively, with a truncation error of 
the order g^^^^- Hence, the expansion parameter in this low-temperature series is T/{gcsjJ,), 
which is also the scaleless parameter appearing in the argument of the leading logarithm 
(remarkably however only after the transverse and the longitudinal contributions have been 
added together). The combination gesf^ is the scale of the Debye mass at high chemical po- 
tential, whose leading-order value is mo = fi'efr/^/Tr. In fact, the calculation of the coefficients 
in (fTT|) required keeping the leading-order "hard-dense-loop" (HDL) part of the gauge boson 
propagator 0,^^, in particular the dynamic screening in but also a HDL correction 
to the real part of the transverse self energy in Q. The above calculation is therefore in 
a certain sense another application of HDL resummation ^^l, which thus turns out to be 
necessary also for a perturbative treatment of the low-temperature regime T//x -C 

As a check on our result and also as a test of its convergence properties, we compare the 
anomalous transverse contributions 5*^,^6 with those of the exactly (albeit only numerically) 
solvable large- A^/ limit [10] in Fig. ^ We find good convergence to the exact result as long 
as T/yU ^ (yfeff/(27r^). This is also the region where the complete laige-Nf result for the 
low-temperature entropy [10] has the anomalous property of exceeding the ideal-gas value. 

Our results do not, however, seem to agree with the results of Ref. 0] which recently 
questioned the presence of a term oc aTlnT"^. The (renormalization group resummed) 
result reported therein rather corresponds to a leading nonanalytic aT^lnT term when 
expanded out perturbatively, which is in fact the type of nonanalytic terms that appear 
already in regular Fermi-liquids 3 • 



For potential phenomenological applications in astrophysical systems, the specific heat 
at constant volume and number density is of more direct interest than the entropy density 
that we have calculated so far. The former (per unit volume) is given by 
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FIG. 2: The perturbative result for the specific heat, normahzed to the ideal-gas value, to order 
2^5/3 g^j^j 2^7/3 ||iQ-^ygp a,nd upper curves, respectively) for two particular values of Os in two- flavor 
QCD (chosen for comparability to Ref. Q) and geflf ~ 0.303 for QED. The deviation of the QED 
result from the ideal-gas value is enlarged by a factor of 20, and the plot terminates where the 
expansion parameter (7r^T)/((/eff/^) ~ 1- 




where Af is the number density, but to the order of accuracy of our expansions, Cy can be 
simply obtained as the logarithmic derivative of the entropy: 



Cv = T —] +OiT'). (15) 




Explicitly, the result is 



+ lc,/3T'^' + lcr/3r/' + 0{T'). (16) 

with C° = NNffj.'^T/3 + 0{T^), and C5/3, C7/3 given by Eqs. dH, (US)). 

For illustrative purposes, we evaluate the ratio of as given by (fTTIjl to the ideal-gas 
value C° for QCD with two massless quark flavors in Fig.|2l using alternatively two values for 
as which have been used also in Ref. It\ and which correspond to one-loop running couplings 
with renormalization point 0.5 GeV (full line) and 1 GeV (dashed line). The shaded bands 
shown are limited from below and above by the results to order T^/^ and T''/'^, respectively, 
and thus indicate the quality of the low-temperature expansion. One may interpret these 
results as roughly corresponding to QCD with a quark chemical potential of 0.5 GeV and the 
total variation corresponding to different renormalization schemes with minimal subtraction 
scale varied between fi and 2/i. The critical temperature for the color superconducting 
phase transition may be anywhere between 6 and 60 MeV so the range T//i > 0.012 in 
Fig. 121 might correspond to normal quark matter. While it is certainly questionable to apply 
perturbative results for as ^ 0.65, Fig. |21 suggests that the anomalous feature of an excess 
of the specific heat over its ideal-gas value may possibly come into pla y in astrophysical 



situations, in particular in the cooling of (proto-)neutron stars [17|, |18|. This should be 
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contrasted with the ordinary perturbative estimate for C^/C^ based on the well-known 11 



exchange term oc g"^ (which, as we have shown, requires T/^^ g). The latter would predict 
CJC^ < 0.6 for a, > 0.65. 

For completeness, we also give the numerical results corresponding to QED, where ges ~ 
0.303. Here the range of temperature, where the specific heat exceeds the ideal-gas value, 
and the deviations from the latter, are much smaller (the deviations from the ideal-gas value 
have been enlarged by a factor of 20 in Fig. El to make them more visible). 

To summarize, we have presented a quantitative evaluation of the leading contributions 
to the entropy and specific heat of high-density QCD and QED in the regime T / ^ <^ g <^ 1, 
which is dominated by non- Fermi-liquid behavior. While the effect remains small in QED, 
it seems conceivable that the anomalous terms in the specific heat play a noticeable role 
in the thermodynamics of a normal quark matter component of neutron or proto-neutron 
stars. 
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